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ABSTRACT 

We present the first high resolution (R=20000-45000, corresponding to 14 
km/s at 4200A to 6.6 km/s at 9000A ) observations of the optical afterglow 
of Gamma Ray Bursts. GRB 020813 and GRB 021004 were observed by 
UVES@VLT 22.19 hours and 13.52 hours after the trigger, respectively. These 
spectra show that the inter-stellar matter of the GRB host galaxies is complex, 
with many components contributing to each main absorption system, and 
spanning a total velocity range of up to about 3000 km/s. Several narrow 
components are resolved down to a width of a few tens of km/s. In the case of 
GRB 021004 we detected both low and high ionization lines. Combined with 
photoionization results obtained with CLOUDY, the ionization parameters of 
the various systems are consistent with a remarkably narrow range with no 
clear trend with system velocity. This can be interpreted as due to density 
fluctuations on top of a regular wind density profile. 



^Based on observations collected at the European Southern Observatory, ESO, the VLT/Kueyen telescope, 
Paranal, Chile, in the framework of programs 69.A-0516(B) and 70.A-0599(B) 
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Subject headings: gamma rays: bursts - cosmology: observations - galaxies: 
abundances - ISM 



1. Introduction 

For a few hours after their onset, Gamma Ray Bursts (GRBs) are the brightest beacons 
in the far Universe, offering a superb opportunity to investigate both GRB physics and high 
redshift galaxies. Tens of minutes after a GRB, its optical afterglow can be as bright as 
magnitude 13-15; a few hours later a magnitude of 16-19 is often achieved. Bright examples 
are the cases of GRB 990123 (z=1.600), for which the reverse shock reached R=9-10 mag 
1 minute after the GRB and the optical afterglow reached R=14 mag 12 minutes after the 
GRB (e.g. Akerlof et al. 1999, Galamaet al. 1999), and GRB 030329 (z=0.1685), for which 
the optical afterglow reached R=12.7 mag at 1.5 hours from the GRB and decreased down 
to R=19 mag after ~ 10 days (e.g. Price et al. 2003; Stanek et al. 2003). High resolution 
(a few tens of km/s in the optical band), high quality (signal to noise > 10 per resolution 
element) spectra can therefore be gathered, provided that the afterglow is observed rapidly 
by 8m class telescopes. Since GRBs are associated with the collapse of massive stars 
(see Woosley 1993, Paczynski 1998, MacFadycn & Wooslcy 1999, Victri & Stella 1998 for 
theoretical reasons and Galama et al. 1998, Stanek et al. 2003, Hjorth et al. 2003 and 
Delia Valle et al. 2003 for observational reasons), this is expected to open a new window in 
the study of the environment in which intermediate to high redshift star-formation occurs 
and, in particular, on the physical, chemical and dynamical state of the inter-stellar matter 
(ISM) of the GRB host galaxies. 

The study of z> 1 galaxies has so far mostly relied upon Lyman-break galaxies (LBGs) 

at z= 3 4 (see e.g. Steidel et al. 1999) and galaxies which happen to be along the hne 

of sight to bright background QSOs. Some of these systems are associated with Damped 
LycK systems (DLA, see e.g. Pettini et al. 1997). However, LBGs are characterized by 
pronounced star-formation and their inferred chemical abundances may be related to these 
regions rather than being representative of typical high z galaxies. Metal line systems 
associated to DLAs along the line of sight to quasars probe mainly galaxy haloes, rather 
than their bulges or discs. Furthermore, it is not clear whether galaxies associated with 
DLAs are truly representative of the whole high-z galaxy population. GRB afterglows 
provide an independent tool to study the ISM of high z galaxies. Savagho, Fall & Fiore 
(2003) studied the metal abundances in three GRB host galaxies using low-medium 
resolution spectroscopy and a curve of growth analysis. They found metal column densities 
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higher that in QSO-DLAs and a strong inverse correlation between [Si/Zn], [Cr/Zn] and 
[Fe/Zn] and the Zn column density, indicating a dense environment and large dust depletion. 
On the other hand, Vreeswijk et al. (2004) found indications for a relatively low metallicity 
and low dust content in the ISM of the z=3.372 host galaxy of GRB 030323, using F0RS2 
low and intermediate resolution spectroscopy. Using again a curve of growth analysis, and 
taking advantage of ultra-deep Gemini multi-object spectrograph observations, Savaglio 
et al. (2004) studied the ISM of a sample of faint K band selected galaxies at 1.4<z<2.0, 
finding Mg II and Fe II abundances much higher than in QSO-DLAs and similar to those 
in GRB host galaxies. A much better job can be done with high resolution {R > 20000) 
observations because: (i) lines can be separated; (ii) metal column densities can be measured 
through a fit to the line profile; (iii) fainter lines can be measured; (iv) information on the 
gas dynamics in the GRB host galaxies can be derived. Furthermore, as suggested by the 
comparison of the Vreeswijk et al. (2004) and the Savaglio et al. (2004) studies, GRB 
afterglows allow us to probe galaxy ISM at much higher redshifts than even ultra-deep, 
standard galaxy spectroscopy, such as the Gemini Deep Deep Survey. 

For all these reasons we started a pilot program to observe bright GRB afterglows of 
promptly localized GRBs with UVES@iVLT. The program has been conceived and designed 
to make full use of the GRB afterglows discovered by the GRB-dedicatcd Swift mission, 
launched on November 20 2004^. Swift will provide GRB positions (< 4' precision) in 
< 10 s. X-ray afterglow positions (3" precision) in < 100 s, and an optical finding chart 
in < 300 s, with a few arcsec to sub-arcsec position accuracy, thus revolutionizing fast 
response multiwavelength studies of GRBs. 

This article concentrates on high resolution spectroscopy of two GRB afterglows. The 
HETE-2 FREGATE, WXM and SXC instruments detected GRB 021004 on 2002 Oct. 4 
12:06:13.57 UT (Shirasaki et al. 2002) and GRB 020813 on 2002 Aug. 13 02:44:19.17 UT 
(Villasenor et al. 2002). The WXM fiight locahzation software produced a rehable and 

relatively accurate position 49 seconds after the burst for GRB 021004 (error box of 30' 
radius) and 4 minutes after the burst for GRB 020813 (error box of 14' radius). Bright 
optical counterparts were identified 10 minutes and 1.9 hours after the triggers by Fox 
et al. (2002, 2003), and Fox, Blake & Price (2002) respectively. Ground analysis of the 
HETE2 data performed a few hours after each GRB improved the localization of the event, 
providing error boxes of 2' and 1' radius, respectively. UVES observations started 13.52 
hours after the GRB 021004 trigger, (Savagho et al. 2002) and 22.19 hours after the GRB 
020813 trigger (Fiore et al. 2002). 



^http://swift.gsfc. nasa . gov 
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Optical low to intermediate resolution spectroscopy of the afterglow of GRB 021004 
was obtained by MoUer et al. (2003), Mirabal ct al. (2003), Matheson et al. (2003), 
Schaefer et al. (2003). One of the main conclusion of these works is that density fluctuations 
on top of a regular wind density profile are able to explain both the presence of strong, 
blue-shifted, high ionization absorption line systems and irregularities in the optical light 
curve. Schaefer et al. (2003) propose a scenario where a clumpy wind originates from a 
massive star progenitor (a Wolf-Rayet star), while Mirabal et al. (2003) suggest that the 
ionized absorption takes place in a fragmented shell nebula around the GRB progenitor. 
Optical spectropolarimetry of the afterglow of GRB 020813 has been presented by Barth et 
al. (2003) who report the detection of strong absorption systems at z=1.223 and z=1.255. 
A further analysis of the Keck LRIS spectrum is presented by Savaglio et al. (2004). The 
present paper reports the results of higher resolution spectroscopy of both afterglows. 

Both GRBs were observed in the X-rays too by the Chandra High Energy Transmission 
Grating System, which has a resolution of ~ 1000 at 1 keV (Butler et al. 2003 and 
references therein). Both observations lasted about 80 ks and started 21 hours and 20 
hours after the GRB event for GRB 020813 and GRB 021004 respectively Fading X-ray 
afterglows were detected in both cases with decay indices of —1.38 ± 0.06 and —0.9 ± 0.1, 
consistent with the values reported for the optical afterglows by Covino et al. (2003) 
and Holland et al. (2003). The time averaged fluxes were 2.2 x 10~^^ erg cm~^ s~^and 
6.3 X 10~^^ erg cm~^ s~^respectively. The X-ray spectra are consistent with power laws of 
energy index of —0.85 ± 0.04 and —1.01 ± 0.08 reduced at low energy by Galactic column 
density (which is rather high in the direction of GRB 020813, NnGai = 7.5 x 10^° cm~^, 
and is Nucai = 4.2 x 10^° cm"^ in the direction of GRB 021004). The upper limits to the 
rest frame X-ray absorbing column are of the order of a few x 10^*^ in both cases. Butler et 
al. (2003) report the detection of an emission line at 1.3 keV at the 3.3cr confidence level, 
interpreted as a SXVI la line blue-shifted by 0.12c. 



2. Observations and data reduction 

In the framework of ESO programs 69.A-0516 and 70.A-0599 we observed the 
afterglows of GRB 020813 and GRB 021004 with the high resolution UV-visual echelle 
spectrograph (UVES, Dckker et al. 2000), mounted at the VLT-UT2 telescope. Table 1 
gives the log of the observations used in this paper. A further 1 hour dichroic 1 observation 
is available in the ESO archive, but we were not able to extract a reliable spectrum from 
this observation, which has therefore been discarded. In order to maximize the signal to 
noise ratio the GGD was rebinned 2x2 pixels. The data reduction has been performed 
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using the UVES pipeline^. The final useful spectra extend from about 4250A to about 
9400A and were rebinned to 0.1 A to increase the signal to noise. The resulting resolution 
element, set to two pixels, ranges then from 14 km/s at 4200A to 6.6 km/s at 9000A . 
The noise spectrum, used to determine the errors on the best fit line parameters, has been 
calculated from the real background subtracted and rebinned spectrum using line free 
regions. This therefore takes into account both statistical errors and systematic errors in 
the pipeline processing and background subtraction. Table 2 gives the signal to noise ratio 
per 0.1 A pixel at different wavelengths for the 4 spectra in Table 1. Figure 1 plots the full 
UVES spectrum of GRB 021004, for which both UVES dichroics, and both red and blue 
arms, were used, allowing us to obtain a particularly wide spectral coverage, extending 
from ~4250A to ^9400 A . A further dichroic 1, blue arm image, covering the band 3400A 
- 3900 A was obtained, but the signal to noise is too low to allow a reliable extraction of 
the spectrum in this band. Figure 2 plots the full UVES spectrum of GRB 020813. This 
spectrum, as well as the spectrum in figure 1, was smoothed with a gaussian function with 
a — 1.5 pixels. Tables 3 and 4 give the equivalent width of all lines detected in the UVES 
spectra for the two GRB afterglows, along with their identification. The equivalent widths 
are computed at the redshift of the absorption systems. For the GRB host galaxy identified 
systems and the 2 main intervening systems at z=1.60 and 1.38 we report also faint lines 
with (signal to noise > 1). Unidentified lines are reported only if the signal to noise is > 3. 
Their equivalent widths are computed at zero redshift. 

3. Column densities 

In this paper we focus on the systems which are likely to be associated with the 
GRB host galaxies. The line fitting was performed using the MIDAS package FITU^MAN 
(Fontana & Ballester 1995). This uses a Voigt profile and yields independently the 
column density N and the Doppler parameter b for each absorption component. For each 



^ http : / / www. eso .org/ observing / dfo /quality /UVES /pipeline/ 



Tabic 1: Journal of observations 



Date UT 


Dichroic 


B. Arm C.W. 


R. Arm C.W. 


slit" 


seeing" 


exposure** 


time since GRB° 


14/08/02 00:55:25 


1 


3460A 


5800A 


1 


<1 


60 


22.19 


05/10/02 01:37:37 


2 


4370A 


8600A 


1 


<1 


30 


13.52 


05/10/02 02:10:53 


1 


3460A 


5800A 


1 


<1 


30 


14.08 


05/10/02 04:09:54 


1 


3460A 


5800A 


1 


<1 


60 


16.06 



"arcsec; ^min; '^hr 




Fig. 1. — UVES spectrum of the GRB 021004 afterglow smoothed with a gaussian function 
with a — 1.5 pixel . The dotted line is the error spectrum 
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Fig. 1. — continued 
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Fig. 1. — continued 
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Fig. 2. — UVES spectrum of the GRB 020813 afterglow smoothed with a gaussian function 
with a — 1.5 pixel . The dotted line is the error spectrum 
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absorption system several lines, spread over the entire spectral range covered by the UVES 
observations, were fitted simultaneously, using the same number of components for each 
line, and the same redshift and b value for each component. 

3.1. GRB 021004 

For this GRB we consider the absorption systems at the following redshifts and 
velocities (in km/s) with respect to the Lya emission of the host galaxy (Mirabal et al. 
2003): z=2.328, v=0 (system A in figure 3); z=2.328, v=-139 km/s (system B); z=2.328, 
v=-224 km/s (system C); z=2.321; z=2.298; z=2.296. For these systems we detected C IV, 
C II, Si IV, Al II, Al III, Mg II and Fe II lines (see Table 3). We report in the table also the 
tentative identification of a Si IIA1304 line at z=2.2953. This is close but not coincident 
with the redshift of systems z=2.296. We therefore consider this identification uncertain, 
and we do not consider this line in the following analysis. 

Unfortunately, the Lja absorption associated with the z=2.328 systems falls exactly 
in the gap between the dichroic 1 and 2 blue arm spectra, and therefore is not covered in 
this study. Zn II and Cr II absorption lines for the main z=2.328 and z=2.321 systems fall 
in the region affected by atmospheric telluric features (see figure 1) and are therefore not 
accessible. 

Figure 3 shows the C IV and Si IV doublets of the above 6 systems. Note that the 
C IVA1548 lines of the z=2.328 A,B,C systems are strongly blended with the C IVA1550 
line of the z=2.321 system, while the Si IVA1404 fine of the z=2.298 system is blended 
with the Si IVA1393 line of the z=2.321 system. Each of the six systems actually comprises 
several components, within a velocity range of several tens of km/s. For this reason the 
identification of the different systems is somewhat subjective, the true message being that 
the geometry and kinematics of the ISM clouds probed by the GRB line of sight are 



Tal)le 2: 


Signal to noise 


ratio 




A 


14/08/02 OOhr 


05/10/02 01hr-02hr 


05/10/02 04hr 


4250A 




3.6 




5000A 


4.7 


5.9 


9.3 


6OOOA 


5.6 


5.4 


11.2 


7OOOA 


5.9 


9.8 


10.2 


7800A 




6.8 




9250A 




2.9 





-n- 



Table 3: GRB 021004 UVES line identifica tions 

A(A) Wresti^T ID Z 



4298, 


.40 


0, 


,13±0, 


,02 


?Si IIA1304.37? 


?2.2953? 


4301, 


.73 


0, 


.06±0, 


.03 


Si IIA1808.01 


1.3793-8 


4304, 


.53 


0, 


.09±0, 


.03 


Si IIA1808.01 


1.3807 


4311, 


.20 


0, 


.16±0, 


.08 






4346 


.30 


0, 


,27±0, 


,08 


Al IIA1670.79 


1.6014 


4347, 


.25 


0, 


.18±0, 


.03 


Al IIA1670.79 


1.6019 


4348, 


.60 


0, 


.11±0, 


.03 


Al IIA1670.79 


1.6028 


4350, 


.95 


0, 


,53±0, 


,08 






4358, 


.28 


0, 


.35±0, 


.08 






4398, 


.34 


0, 


.06±0, 


.03 


C IIA1334.53 


2.2958 


4401, 


.68 


0, 


.24±0, 


.02 


C IIA1334.53 


2.2981-4 + 












?C IPA1335.70? 


?2.2958? 


4415, 


.50 


0, 


.10±0, 


.03 


Al IIIA1854.71 


1.3807 


4441, 


.34 


0, 


.04±0, 


.02 


C IIA1334.53 


2.328 


4445, 


.35 


0, 


.09±0, 


.02 


C IPA1335.70 


2.328 


4593, 


.59 


0, 


,15±0, 


,02 


Si IVA1393.76 


2.2958 


4596, 


.71 


0, 


.02±0, 


.02 


Si IVA1393.76 


2.2981-4 


4623, 


.48 


0, 


.21±0, 


.02 


Si IVA1402.77 


2.2958 


4626, 


.41 


0, 


.02±0, 


.02 


Si IVA1402.77 


2.2981-4 


4629, 


.72 


0, 


,40±0, 


,02 


Si IVA1393.76 


2.321 


4631, 


.35 


0, 


.30±0, 


.07 






4635, 


.85 


0, 


.65±0, 


.07 






4637, 


.38 


3, 


,63±0, 


,07 






4638, 


.01 


0, 


,31±0, 


,02 


Si IVA1393.76 


2.328 


4643, 


.19 


0, 


.25±0, 


.07 






4658, 


.07 


0, 


.27±0, 


.02 


Si IVA1402.77 


2.321 


4666, 


.28 


1, 


,07±0, 


,07 






4668, 


.38 


0, 


.56±0, 


.02 


Si IVA1393.76 


2.328 


5102, 


.54 


0, 


.50±0, 


.02 


C IVA1548.20 


2.2958 


5106, 


.50 


0, 


.35±0, 


.02 


C IVA1548.20 


2.2981-4 


5111, 


.04 


0, 


,33±0, 


,02 


C IVA1550.77 


2.2958 


5114, 


.82 


0, 


.16±0, 


.02 


C IVA1550.77 


2.2981-4 


5141, 


.11 


1, 


.74±0, 


.02 


C IVA1448.20 


2.321 


5150, 


.4 


2, 


.10±0, 


.02 


C IVA1550.77 


2.321 + 












C IVA1448.20 


2.328 


5159, 


.66 


1, 


.57±0, 


.02 


C IVA1550.77 


2.328 


5510, 


.62 


0, 


.05±0, 


.02 


Al IIA1670.79 


2.2981-4 


5560 


.30 


0, 


,47±0, 


,02 


Al IIA1670.79 


2.328 


5577, 


.71 


1, 


.47±0, 


.03 


Fe IIA2344.21 


1.3793-8 


5580, 


.81 


0, 


.07±0, 


.03 


Fe IIA2344.21 


1.3807 


5592, 


.94 


0, 


.14±0, 


.06 






5649 


.59 


0, 


,31±0, 


,03 


Fc IIA2374.21 


1.3793-8 


5669, 


.68 


0, 


.67±0, 


.03 


Fe IIA2382.77 


1.3793-8 


6064, 


.56 


0, 


.60±0, 


.06 






6098, 


.32 


0, 


,22±0, 


,06 







Errors are 67% confidence intervals. the equivalent width of unidentified lines is computed 
at zero redshift; ? = uncertain identification. In many cases several components contribute 
to a single line in one entry of the table. In these cases either the redshift is given with only 
4 Hiprits nr it i.s p-ivftn a.t; a, ra,np"ft A -I- sipTi inHira.tftt; tha.t a, linp is nartiv biftnrlftrl with thft 



- 12 - 



Table 3: GRB 021004 UVES line identifications, continued. 



ApO Wreet{k )« ID 



6099, 


.33 


0, 


,19±0, 


,02 


Fc IIA2344.21 


1.6019 


6101, 


.32 


0, 


,09±0, 


,02 


Fc IIA2344.21 


1.6028 


6154, 


.64 


0, 


.52±0, 


.03 


Fe IIA2586.65 


1.3793-8 


6158, 


.03 


0, 


.01±0, 


.03 


Fe IIA2586.65 


1.3807 


6160 


.76 


0, 


,17±0, 


,06 






6172, 


.42 


0, 


,35±0, 


,02 


Al IIIA1854.72 


2.328 


6176, 


.74 


0, 


.21±0, 


.06 






6177, 


.88 


0, 


.17±0, 


.02 


Fe IIA2374.46 


1.6019 


6180, 


.04 


0, 


,14±0, 


,05 






6186, 


.51 


0, 


.25±0, 


.02 


Al IIIA1862.79 


2.321 + 


6186, 


.99 


0, 


.11±0, 


.02 


FeIIA2600.17 


1.3793-8 


6187, 


.81 


0, 


.05±0, 


.02 


FeIIA2600.17 


1.3798 


6190 


.23 


0, 


,12±0, 


,02 


Fe IIA2600.17 


1.3807 


6199, 


.36 


0, 


.11±0, 


.02 


Al IIIA1862.79 


2.328 + 


6198, 


.4 












6199, 


.53 


0, 


,29±0, 


,02 


Fc IIA2382.77 


1.6019 


6201, 


.66 


0, 


,28±0, 


,02 


Fc IIA2382.77 


1.6028 


6653, 


.88 


0, 


.89±0, 


.02 


Mg IIA2796.35 


1.3793-8 


6655, 


.88 


0, 


.24±0, 


.02 


Mg IIA2796.35 


1.3802 


6657, 


.60 


0, 


,47±0, 


,02 


Mg IIA2796.35 


1.3807 


6671, 


.00 


0, 


.78±0, 


.02 


Mg IIA2803.53 


1.3793-8 


6672, 


.73 


0, 


.06±0, 


.02 


Mg IIA2803.53 


1.3802 


6674, 


.77 


0, 


.32±0, 


.02 


Mg IIA2803.53 


1.3807 


6730 


.16 


0, 


,25±0, 


,02 


Fc IIA2586.65 


1.6019 


6732, 


.55 


0, 


.09±0, 


.02 


Fe IIA2586.65 


1.6028 


6750, 


.64 


0, 


.08±0, 


.02 


Mn IIA2594.50 


1.6019 


6764 


.18 


0, 


,43±0, 


,04 






6765, 


.43 


0, 


,25±0, 


,02 


Fc IIA2600.17 


1.6019 


6767, 


.79 


0, 


.17±0, 


.02 


Fe IIA2600.17 


1.6028 


6768, 


.78 


0, 


.18±0, 


.04 






6781, 


.81 


0, 


,04±0, 


,02 


Mn IIA2606.46 


1.6019 


6788, 


.06 


0, 


.03±0, 


.02 


Mg IA2852.94 


1.3793-8 


7003, 


.91 


0, 


.12±0, 


.03 






7041, 


.91 


0, 


.12±0, 


.03 






7042, 


.52 


0, 


.11±0, 


,03 






7043, 


.25 


0, 


.10±0, 


.03 






7274, 


.52 


0, 


.26±0, 


.01 


Mg IIA2796.35 


1.6014 


7275, 


.83 


0, 


.33±0, 


.01 


Mg IIA2796.35 


1.6019 


7277, 


.28 


0, 


,49±0, 


,01 


Mg IIA2796.35 


1.6024 


7278, 


.39 


0, 


.24±0, 


.01 


Mg IIA2796.35 


1.6028 


7293, 


.12 


0, 


.27±0, 


.01 


Mg IIA2803.35 


1.6014 


7294, 


.58 


0, 


,36±0, 


,01 


Mg IIA2803.53 


1.6019 


7295, 


.91 


0, 


,04±0, 


,01 


Mg IIA2803.53 


1.6024 


7296, 


.85 


0, 


.23±0, 


.01 


Mg IIA2803.53 


1.6028 


7422, 


.05 


0, 


.04±0, 


.01 


Mg IA2852.94 


1.6014 


7423, 


.01 


0, 


,23±0, 


,01 


Mg IA2852.94 


1.6019 



" the equivalent width of unidentified lines is computed at zero redshift. ? = uncertain 
identification. In many cases several components contribute to a single line in one entry of 
the table. In these cases either the redshift is eiven with onlv 4 dieits or it is eiven as a 
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Table 3: GRB 021004 UVES line identifications, continued. 



A(A) 




ID 


z 


7929.74 


0.75±0.05 






7929.53 


0.08±0.02 


Fc IIA2382.77 


2.328 


7930.42 


0.04±0.02 


Fe IIA2382.77 


2.328 


7931.03 


0.03±0.02 


Fe IIA2382.77 


2.328 


8228.16 


0.37±0.05 






9178.45 


0.75±0.08 






9205.45 


0.70±0.09 






9216.41 


0.04±0.03 


Mg IIA2796.35 


2.2958 


9222.70 


0.09±0.03 


Mg IIA2796.35 


2.2981 


9223.62 


0.18±0.03 


Mg IIA2796.35 


2.2984 


9246.33 


0.09±0.03 


Mg IIA2803.53 


2.2981 


9247.39 


0.12±0.03 


Mg IIA2803.53 


2.2984 


9330.33 


0.57±0.03 


Mg IIA2803.53 


2.328 


" the equivalent width of unidentified lines is computed at zero redshift 



complex. Nevertheless, sticking to the above system identifications will be useful. 

There are 13 lines in Table 3 associated with the z=2.328 systems, 5 lines associated 
with the z=2.321 system, and up to 17 lines associated to the z=2. 296-2. 298 systems. 
Some of these lines can be split further in several components. To test the robustness of 
the fit, in terms of accuracy and stability of the results, we performed many fits, using 
several combination of lines/systems. We found that the fits presented below are a good 
compromise between increasing the statistical precision of the fit, obtained increasing the 
number of lines fitted simultaneously, and the stability/repeatibility of the results, which 
degraded increasing the number of fitted parameters, due to the increasingly complex shape 
of the hyper-surface in the parameter space, that may contain many local minima. 

For the z=2.328 A,B,C and z=2.321 systems we fitted simultaneously the C 1VAA1550, 
1548, C 11A1334, C 11*A1335, Si IVAA1404, 1393 Al 11 A1670 and the Al IIIA1854 line. We 
excluded from the fit the z=2.328 Al IIIA1862 line because its blue wing is strongly blended 
with another line, and the z=2.321 Al III A1862 line because it is strongly blended with 
the Fe IIA2600 line of one of the components of the z=1.38 intervening system. For the 
same systems we also fitted the Fc I1A2382, Fc IIA2374 Fc IIA2344 Fc IIA1608, Mg I1A2803 
and Mg IIA2796 lines simultaneously. Figure 4 shows the line spectra in velocity space, 
along with the best fit model, while Table 5 presents the best fit abundances along with 
the velocity shift of each system with respect to the redshift of the host galaxy, assumed to 
be 2.328 for GRB 021004 (Mirabal et al. 2003). We used 2 components for the z=2.328_A 
system, and 1 component for the other 3 systems. The best fit doppler parameter b ranges 
from a minimum of 12±6 km/s (one of the components of the z=2.328_A system) up to 
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Table 4: GRB 020813 UVES line identifications 



A(A) 


Wrest (A 




ID 


z 


4824, 


.64 


1, 


,06± 


0, 


,05 






4824, 


.64 


0, 


.59± 


0, 


,05 






4998, 


.03 


0, 


.47± 


0. 


.05 






5212, 


.23 


0, 


.18± 


0. 


.03 


Fe IIA2344.21 


1.2234 


5236, 


.74 


1, 


,64± 


0, 


,07 






5281, 


.33 


0, 


.09± 


0, 


,03 


Fe IIA2344.21 


1.2529 


5285, 


.19 


1, 


.70± 


0, 


,03 


Fe IIA2344.21 


1.255 


5297, 


.11 


0, 


,85± 


0, 


,03 


Fc IIA2382.77 


1.2234 


5353, 


.61 


1, 


.31± 


0, 


,03 


Fe IIA2374.46 


1.255 


5368, 


.21 


0, 


.28± 


0, 


,03 


Fe IIA2382.77 


1.2529 


5372, 


.23 


1, 


.63± 


0. 


.03 


Fe IIA2382.77 


1.255 


5840, 


.73 


0, 


,44± 


0, 


,06 






5850, 


.86 


0, 


.73± 


0. 


.06 






5858, 


.24 


0, 


.23± 


0. 


.02 


Fe IIA2600.17 


1.2529 


5862, 


.33 


1, 


.82± 


0. 


.02 


Fe IIA2600.17 


1.255 


5896, 


.25 


0, 


,40± 


0, 


,06 






6067, 


.05 


2, 


.30± 


0, 


,06 






6216, 


.93 


1, 


.34± 


0. 


.03 


Mg IIA2796.35 


1.2234 


6232, 


.85 


1, 


.36± 


0, 


,03 


Mg IIA2803.53 


1.2234 


6300, 


.13 


0, 


,55± 


0, 


,03 


Mg IIA2796.35 


1.2529 


6304, 


.72 


2, 


.38± 


0, 


,03 


Mg IIA2796.35 


1.255 


6310, 


.38 


0, 


.30± 


0, 


,06 






6316, 


.13 


0, 


,30± 


0, 


,03 


Mg IIA2803.53 


1.2529 


6321, 


.11 


2, 


,19± 


0, 


,03 


Mg IIA2803.53 


1.255 


6343, 


.49 


0, 


.11± 


0, 


,03 


Mg IA2852.13 


1.2234 


6427, 


.36 


0, 


.24± 


0. 


.03 


Mg IA2852.13 


1.2529 


6432, 


.12 


1, 


,34± 


0, 


,03 


Mg IA2852.13 


1.255 


6686, 


.87 


0, 


,27± 


0, 


,03 


Fc IA2967.76 


1.2529 



Errors are 67% confidence intervals. " the equivalent width of unidentified lines is computed 
at zero redshift. 

a maximum of 107±10 km/s (system z=2.321). Of course different b values are obtained 
using an higher (or lower) number of components. Unfortunately the signal to noise ratio of 
our spectrum is not good enough to unambiguously identify the different components. On 
the other hand, we verified that the total best fit column density of each system is stable, 
within the statistical errors, changing the number of components in each system. 

Similar series of fits were performed for the z=2.296 and z=2.298 systems, for which 
we used 1 and 2 components, respectively (sec figure 5). Table 5 gives again the best fit 
abundances for these 2 systems. The best fit doppler parameter b ranges from 9±7 km/s 
to 40±13 km/s. 
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Fig. 3. — C IV (upper panel) and Si IV (lower panel) absorption systems in the UVES 
spectra of GRB 021004 for z=2.328 (see Table 5). 
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Fig. 4.— GRB 021004. Upper panel: UVES spectrum near the C IVAA1550, 1548, C IIA1334, 
C ir A1335, Si IVAA1404, 1393, Al IIA1670 and Al IIIA1854 lines for the z=2.328 and z=2.321 
systems in velocity space, along with the best fit model (solid line) and residuals (grey). 
Lower panel: same for the Fe IIA2382, Fe IIA2374, Fe IIA2344, Fe IIA1608 and Mg IIA2803 
lines. The zero in the velocity scale refers to the redshift of the host galaxy. 
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Fig. 5.— GRB 021004. Upper panel: UVES spectrum near the C IV AA 1550, 1548, 
C IIA1334, C II*A1335, Si IVAA1404, 1393 lines for the z=2.298 and z=2.296 systems in 
velocity space, along with the best fit model (solid line) and residuals (grey). Lower panel: 
same for the Fe IIA2382, Fe IIA2374, Fe IIA2344, Fe IIA1608, Mg IIA2803 and Mg IIA2796 
lines. The zero in the velocity scale refers to z=2.296. 
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The detection of both high and low ionization hnes, feasible thanks to the extremely 
wide spectral coverage achieved by UVES, allows us to obtain constraints on the ionization 
status of the gas responsible for the UV absorption, by comparing ion column density ratios 
with the predictions of photoionization codes. Unfortunately we were not able to measure 
column densities of different ions of the same element, except for C II/C IV and Al II/ Al III 
for system z=2.328_A, and C II/C IV for system z=2.298. In the other cases we are forced 
to use ratios of column densities of different ions of different elements in this analysis. 
These estimates of the gas ionization parameter are therefore somewhat degenerate with 
respect to relative element abundances. We used the Grevesse & Anders (1989) meteoritic 
abundances with extensions by Grevesse et al. (1993). 

Figure 6 shows the logarithmic ratio between the C IV, C II, Fe II, Al II and Mg II 
column densities and that of Si IV for the six absorption systems as a function of the 
velocity shift with respect to the redshift of the host galaxy. No large variation of the 
ion ratios are seen for the six systems. We compared these line ratios to the predictions 
obtained by simulating a gas cloud illuminated by an ionizing continuum. We used Cloudy 
(vs 90.04, Ferland et al. 2002) to build grids of photoionization models as a function of U, 
the ionization parameter. U is defined as the ratio between the ionizing photon density 
and the electron density of the gas. It is computed assuming a constant density profile 
throughout the cloud, and a plane parallel geometry. We studied gas densities between 1 
cm~^ and 10^ cm~^. The ionizing continuum was assumed to be a power law, F{E) — E"^ 
photons cm~^ s~^, with cutoffs at low and high energies. The high energy cutoff is fixed at 
lO^^Hz, while we run the simulations for different low energy cutoffs, from 10^°Hz to lO^^Hz. 
We produced several sets of simulations with F in the range 1-2. The ionizing continuum 
is constant in time. GRBs are highly variable sources and the ionization structure of the 
gas can be better studied using a time-dependent photoionization code, such as those of 
Nicastro et al. (1999) and Perna & Lazzati (2002). Nevertheless, our simpler approach is 



Table 5: GRB 021004. Logarithmic ion column densities in cm ^ 



System 


V. shift'' 


Si IV 


CIV 


CII 


CII* 


Al II 


Fell 


Mgll 


2.328_A 





15.30±0.56 


>15.2 


13.40±0.45 


13.90±0.58 


13.55±0.35'= 


13.34±0.15 


13.82±0.24 


2.328_B 


-139 


14.27±0.16 


>14.4 


<13.2 


<13.2 


12.23±0.35 


12.53±0.24 


12.93±0.60 


2.328_C 


-224 


13.24±0.17 


14.40±0.11 


<13.2 


<13.2 


<12.0 


<12.3 


<12.8 


2.321 


-632 


14.11±0.08 


15.04±0.05 


<13.2 


<13.2 


<12.0 


<12.3 


<12.8 


2.298 


2729 


13.43±0.33 


14.16±0.20 


14.28±0.20* 


<13.2 


<12.0 


12.68±0.20 


13.22±0.10 


2.29(i 


2913 


13.79±0.1() 


14.72±0.2o 


13.43±U.40 




<12.() 


<12,3 


12.()8±0.21 



errors, upper and lower hmits are 90% confidence intervals; "km/s; the z=2.296 C 11* fine 
is completely blended with the z=2.298 C II line; ^ for system z=2.328_A we measuered also 
a column density of Al III of 13.66±0.30. 
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Fig. 6. — Tlie logaritlimic ratio between tlie C IV, C II, Fe II, Al II and Mg II column 
densities to tliat of Si IV for the six absorption systems as a function of the negative velocity 
shift with respect to the redshift of the host galaxy. 



instructive in identifying general trends. 

The predicted photoionization equihbrium ion ratios of [Mg Il/Si IV], [Fe Il/Si IV], 
[Al Il/Si IV] and [C IV/Si IV] are shown in figure 7 for F = 2 and n=l cm~^ and figure 8 for 
F = 2 and n=10^ cm~^. For each of the six redshift systems, thick segments superimposed 
on these curves indicate our 90% confidence determinations for these ratios. In figure 7 
(for n=l cm-3) the [Mg Il/Si IV], [Fe Il/Si IV] and [Al Il/Si IV] ion ratios for the systems 
z=2.328_A, z=2.328_B, z=2.328_C, and z=2.296 are all consistent with the same ionization 
status. For the system z=2.328_A we measured the column density of Al II and Al III. Also 
their ratio is consistent with the ionization parameter inferred from the above ion ratios, 
see figures 7 and 8. For the same system we measured the column density of C II. The 
[C Il/Si IV] ion ratio indicates a ionization parameter somewhat higher than that indicated 
by the other ion ratios for n=l cm~^ (see figure 7), while for n=10® cm~^ the ionization 
parameter is fully consistent with that indicated by the other ion ratios (figure 8). The 
system at z=2.298 shows a ionization parameter lower than that of system z=2.328_A, 
while the system at z=2.321 shows a higher ionization parameter. However, all ionization 
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Fig. 7. — Ion ratios from CLOUDY models for the six absorption systems, assuming a gas 
density of 1 cm~^: [Mg Il/Si IV] = solid red curves, [Fe Il/Si IV]= dashed blue curves, 
[Al Il/Si IV] = dotted magenta curves and [C IV/Si IV] and [C Il/Si IV]= dot-dashed black 
curves. The line indicating the model values is thick when it is consistent within 90% with 
the measured ion ratios. The vertical solid lines mark the range of allowed U values in each 
panel. 
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Fig. 8. — Ion ratios from CLOUDY photoionization models for the six absorption systems 
assuming a gas density of 10^ cm~^. Curves as in figure 7. The hne indicating the model 
values is thick when it is consistent within 90% with the measured ion ratios, as in figure 7. 
The vertical solid lines mark the allowed U values in each panel. 
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Fig. 9. — The logarithm of the best fit ionization parameters from figure 7 (n=l cm~^ filled 
circles) and from figure 8 (n=10^ cm~^, open triangles, shifted by logU=-1.5), as a function 
of the velocity shift with respect to the redshift of the host galaxy. Error-bars and upper 
limits represent the purely statistical 90% confidence intervals. 
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parameters systems are consistent with a remarkably narrow range: 10~^'^ < U < 10~^, for 
n=l cm~^ and 10^^'^ < U < 10~^, for n=10^ cm~^ see figure 9, which plots the ionization 
parameters as a function of the velocity for the six systems. We note that in these ranges of 
U the [Fe II/Mg II] ratio varies by less than 30%, and the value found for the three systems 
in which we detect Fe II and Mg II lines ([Mg Il/Fe II] f» 0.4) is consistent with what is 
expected from meteoritic abundances (and as was adopted for the CLOUDY calculations). 
This value is also in the range found for the Galactic ISM 0.1<[Mg/Fe]<0.8. It is somewhat 
higher than that found by Savaglio et al. (2004) in a sample of faint K-band selected 
galaxies at 1.4<z<2.0 -0.84<[Mg/Fe]<0.13. Assuming meteoritic abundances the hydrogen 
column density of the densest system (z=2.328_A) is ^5 x 10^^ cm~^, where the equality 
yields if all Si is in the form of Si IV. This column is so low that would have easily escaped 
detection in the Chandra X-ray spectrum (Butler et al. 2003) and it is consistent with the 
upper limit given by MoUer et al. (2003). 

Qualitatively similar results are drawn from the analysis of the n=10^ cm~^ curves, 
although the range of acceptable ionization parameters for the 6 systems is narrower, see 
figure 8. Qualitatevely similar results were also obtained for the F = 1 models. 

3.1.1. The z=2. 328_A system 

For this system we detect the C irA1335 fine structure hue. Unfortunately, the 
C IIA1334 ground state line is barely visible (with a signal to noise of just 2). The presence 
of a strong C 11* line suggests either a high gas density or a strong radiation field (Srianand 
& Pctitjean 2000, 2001, Silva & Viegas 2002). C 11* and Si II* fine structure lines have 
been detected in two other GRBs: GRB 020813 (Savagho & Fall 2004) and GRB 030323 
(Vreeswijk et al. 2004). Their ratio can be used to constrain the gas density or the radiation 
field (if both lines are detected one may solve for both variables, as shown by Srianand & 
Petitjean 2000), so their detection in GRB spectra looks very promising to constrain the 
physical properties of the absorbing gas. 

For this system the C IV line is strongly saturated and so it cannot be used together 
with the C II to constrain the gas ionization status. As discussed in the previous section, 
within rather large uncertainties, the [G Il/Si II] would suggest for n=l cm~^ a ionization 
parameter somewhat higher than that indicated by all other line ratios, assuming meteoritic 
abundances. This may be explained by a slight under- abundance of G with respect to Si in 
this system. A similar conclusion applies to system z=2.321, see below. 

The G IVA1548 fines of the z=2.328 A,B and G systems are blended with the G IVA1550 
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of the z=2.321 system, suggesting the presence of hne locking, as first reported by Savagho 
et al. (2002). Line locking is usually interpreted as the signature of line radiation pressure 
acceleration (see e.g. Foltz et al. 1987, Srianand & Petitjean 2000). 

3.1.2. The z=2.321 system 

For this system (as for z=2.328_C) we were able to detect C IV and Si IV absorption 
hues only. The upper hmits on the Mg II and Fe II transition allow us to put a significant 
lower limit on the ionization parameter of UJ^IO"^ for this system for a gas density of 1 
cm~^ and 10~^'^ for a gas density of 10^ cm^^. For n=l, this is higher than the values for 
the other systems at a confidence level better than 90%. The [C IV/Si IV] ratio is formally 
inconsistent with this higher ionization parameter, but the disagreement is marginal, taking 
into account that the C IV/Si IV curves in figures 7 and 8 are very fiat, and therefore 
small differences in the obsterved C IV/Si IV ratio would translate in large differences in 
the ionization parameter. Also in this case the disagreement may be explained by a slight 
under-abundance of C with respect to Si. 

3.1.3. The z^2. 296 and z^2. 298 systems 

Given the relatively high velocity shift of about 3000 km/s between the z=2.296 and 
z=2.298 systems and the host galaxy redshift, there is the possibihty, at least in principle, 
that these systems are not associated with the GRB host galaxy but rather are intervening 
systems. We evaluated the probability of finding intervening C IV systems within 3000 
km/s from the GRB host galaxy, from the density distribution of C IV systems given in 
D'Odorico et al. (1998) and Petitjean & Bergeron (1994). The probability of finding by 
chance two systems with column densities equal or larger than those given in Table 5 for 
these two lowest redshift systems is 1.2 x 10^^, while the probability for only one of these 
systems is 0.1% for the z=2.296 system and 0.5% for the z=2.298 system. 

The z=2.298 system has a ionization parameter < 10"^ '' at the 90% confidence level 
for a gas density of 1 cm^'^ (around 10~^'^ for a gas density of 10® cm"^). This is the lowest 
of the ionization parameters found in the 6 systems. Furthermore, for this system we detect 
a rather strong C IIA1334 ground state line but not the C irA1335 high excitation line, see 
figure 5, again indicating a lower radiation field (or a lower gas density). 

The observed wavelength of the z=2.298 system C IIA1334 line coincides with the 
wavelength of the z=2.296 C IFA1335 line (see Table 3). This may be a further indication 
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of line locking. 



3.1.4- Intervening systems 

Two main intervening systems are present along the line of sight to GRB 021004, one 
around z=1.60 and the other around z=1.38 (see eg. Mirabal et al. 2003, MoUer et al 
2003). The UVES high resolution spectrum of GRB 021004 allows an accurate study of 
these systems, which are spht in several components (see Table 3). 

The system at z=1.38 consists at least 3 components at z=1.3807 (z=1.38_A), z=1.3802 
(z=1.38_B), z=1.3795 (z=1.38_C). The latter system is also split in several parts in the range 
1.3793-1.3798. The total velocity range spanned by the system is ~ 180 km/s. We detected 
the following hues belonging to this system: Si IIA1808, Al IIIA1854, Fe IlA2344,Fe IIA2374, 
Fe IlA2586,Mg IIAA2796,2803, and Mg IA2852. 

The system at z=1.60 is made by at least 4 components at z=1.6028 (z=1.60_A), 
z=1.6024 (z=1.60_B), z=1.6019 (z=1.60_C) and z=1.6014 (z=1.60_D). The total velocity 

range spanned by the system is of ~ 160 km/s. For this system we detected the 
following Uncs: Al IIA1670, Fc IlA2344,Fe IIA2374, Fc IIA2382, Fe IlA2586,Mn IIA2594, 
Fe IlA2600,Mn IIA2606, Mg IIAA2796,2803, and Mg IA2852. 

We estimated the column densities of the ions of the intervening systems using the 
same fitting procedure used for the six systems associated with the GRB host galaxy. They 
are given in Table 6 

Table 6: Log. ion column densities in cm~^ of the GRB 021004 intervening systems 



System 


V. shift" 


Fell 


Mg II 


Si II 


Al III 


z=1.38_A 





15.68±0.16 


13.23±0.07 


15.12±0.20 


13.02±0.20 


z=1.38_B 


-63 




12.92±0.32 






z=1.38_C 


-150 


14.75±0.40 


15.20±0.70 






System 


V. shift" 


Fell 


Mgl 


Mg II 


Al II 


z=1.60_A 





13.69±0.16 


11.33±0.26 


13.54±0.21 


12.66±0.27 


z=1.60_B 


-46 


12.63±0.16 


11.64±0.18 


12.95±0.10 




z=1.60_C 


-104 


14.35±0.18 


12.80±0.27 


15.45±0.80 


13.24±0.50 


z=1.60_D 


-160 


13.57±0.10 


11.82±0.14 


13.47±0.06 


12.58±0.25 



errors, upper and lower limits are 90% confidence intervals; "km/s 
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3.2. GRB 020813 

For this GRB we considered the absorption systems at the following velocities with 
respect to the redshift of the host galaxy, which was assumed to be 1.255 (Barth et al. 
2003): v=0 km/s (system A in figure 10), v=-106 km/s (system B) and v=-306 km/s 
(system C). Note that the system A and system B lines are strongly blended with the 
exception of the Fe II2374 hne. We also consider the system at z=1.2234. 

We fitted simultaneously the 6 lines of the 3 absorption systems in figure 10. Table 7 
gives for the 3 systems the best fit Fe II and Mg II abundances along with the velocity shift 
of each system with respect to the redshift of the host galaxy. Mg II hnes of systems A and 
B are strongly saturated and therefore their column density estimates are more uncertain. 
Unfortunately the UVES spectrum covers a wavelength range much smaller than the Keck 
LRIS spectrum and several of the lines studied by Barth et al. (2003) and Savaglio & Fall 
(2004), in particular Zn II, Cr II, Si II, are not accessible. 

We performed similar fits to the 5 hnes associated with the z=1.2234 system, spht in 
two components. The results are in figure 11 and in Table 7. 

Unfortunately in this case the redshift of the GRB is not high enough to have strong 
high ionization lines in the spectrum, and therefore we do not have a direct way to constrain 
the ionization status of the gas responsible for the UV absorption. 

The [Fe II/Mg II] ratio is consistent with a constant value in the four systems. This 
ratio is also consistent with meteoritic abundances, assuming that the ionization parameter 
is in the range 10~^-^ — 10~^ 

3.3. The z=1.2234 system 

The system at z= 1.2234 is shifted by about 4200 km/s from the redshift of the host 
galaxy. The probability to find an Mg II intervening system with Wa2796.35 ~ 1-3 ± 0.03A 

Table 7: GRB 020813, Logarithmic ion column densities in cm~^ 



System 


velocity shift" 


Fell 


Mg II 


z=1.255_A 


0. 


15.20±0.29 


> 16.0 


z=1.255_B 


106 


15.47±0.27 


> 15.5 


z=1.255_C 


306 


13.58±0.10 


13.79±0.21 


z=1.2234 


-4204 


13.96±0.18 


> 16.0 



"km/s 
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(see Table 7) within this velocity range is ^1%, using the distribution of Mg II systems 
given by Steidel & Sargent (1992). Although not as conclusive as in the case of the z=2.296 
and z=2.298 systems in the spectrum of GRB 021004, we remark that this probability is 
rather small, and that the velocity shift with respect to the redshift of the host galaxy is 
intriguingly similar to that of the GRB 021004 systems. 



4. Discussion 

GRB 021004 shows absorption systems which span ~ 3000 km/s in velocity towards 
the observer (no receding system is detected). The systems are most likely local to the 
GRB since the probability of finding two or more different absorption systems along the 
line of sight by random fluctuations is negligible. 

It has been suggested that radiative acceleration by the prompt GRB emission may be 
responsible for the high detected speeds (Schaefcr ct al. 2003). However, this seems unlikely 
since low ionization ions such Fe II and Mg II arc present in the high velocity absorber. 
The radiative acceleration is dominated by the GRB and early afterglow radiation. It 
is therefore unlikely that recombination can be so rapid to have any influence on the 
acceleration of the absorbing material. In the photoionization phase each nucleon receives 
an amount of energy comparable to the ionization potential of the X-shell electron. For a 
solar metallicity plasma the acceleration is regulated by H atoms and, at the end of the 
photoionization, the absorbing material has an outward velocity of: 



\ 



'50 km/s (1) 



nip 



Once the gas is fully ionized, acceleration proceeds due to radiation pressure on free 
electrons through Inverse Compton (IC) interactions. For a burst of isotropic equivalent 
energy Eiso the IC acceleration is obtained by momentum conservation: 

vic = /;f^ - 0.6 R^i km/s (2) 
Att nip c 

we conclude that at the minimum distance of the absorbing medium it! ~ 10^^ cm (Lazzati 
et al. 2002; Heyl & Perna 2003; in order to be ahead of the fireball at the time at 
which absorption is detected) radiative acceleration is unable to propel the absorber to an 
outfiowing speed comparable to the value(s) observed. 

Possible explanations for the large outfiowing speed are either a supernova exploded 
several years prior to the GRB (Vietri & Stella 1998) or a high velocity wind from the 
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Fig. 10.— UVES spectrum of GRB 020813 near the Fe IIA2600, Fe IIA2382, Fe IIA2374, 
Fe IIA2344, Mg IIA2803 and Mg IIA2796 lines for the three absorption systems, in velocity 
space, along with the best fit model, solid line and residuals. The zero of the velocity scale 
refers to z= 1.2545 
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Fig. 11.— UVES spectrum of GRB 020813 near the Fe IIA2382, Fe IIA2344, Mg IIA2803, 
Mg IIA2796 and Mg IA2852 lines for the z=1.2234 absorption systems, in velocity space, 
along with the best fit model, solid line and residuals. The zero of the velocity scale refers 
to z=1.2234 
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progenitor Wolf-Rayet (WR) star (MacFadyen & Woosley 1999; Schaefer et al. 2003; 
Mirabal et al. 2003). The large velocity spread with similar ionization parameter is however 
difficult to account for in a SNR scenario, where values around a typical expansion velocity 
would be expected. On the other hand, WR winds are known to be clumpy and velocities 
up to ~ 4000 km/s were detected from P-Cygni proffies (Niedzlieski & Skorzyhski 2002). 

In this work we derived physical parameters of the absorbers under the assumption 
of equilibrium conditions which, ~ 0.5 days after the GRB explosion, are attained only if 
the electron density in the absorber is n ~ 10^^*^ cm^^. Assuming equilibrium conditions, 
the total column density of the absorber can be computed from the ion column density 
corrected for the ionization fraction. Consider the Si IV line in the A system. The density 
is given by: 

]\j 1 r|15.3 

n = ^5x10' R7} ri'l cm~^ (3) 

AsiRf] 4 X lO-noi^i^iglO-Ve 

where R is the fireball radius at the time the lines were observed, Agi is the Si abundance 
and 1] is the ratio of the width of the absorbing shell over its radius. The ionization 
parameters derived above and implications discussed here are therefore relevant only in the 
case of an extremely clumpy wind or SNR. 

The photoionization results of CLOUDY yield an ionization parameter constrained in 
a relatively small range 10~^'^ < U < 10~^. In a single explosion GRB model, the ionization 
parameter scales with the square of the outflow velocity, since the absorber's distance 
scales with velocity as well. Even though the ionization parameter depends on the electron 
density, it seems unlikely that density variations are such as to compensate for the large 
velocity difference. In a wind environment, on the other hand, the ionization parameter 
is constant since the photon density and particle density scale with the same power of the 
distance. The relatively small variations in the inferred U, which do not show any clear 
trend with velocity, can therefore be interpreted as density fluctuations on top of a regular 
R^^ wind density proffie (as already discussed by e.g. Schaefer et al. 2003 and Mirabal et 
al. 2003). 

Finally, we note that the Fe II and Mg II column densities found in GRB 020813 
are 10-100 times higher than those in GRB 021004. This is likely due to a much higher 
ionization of the gas in the latter case, rather than a large difference in the total absorbing 
column of gas or to highly non solar metal abundances (note that in both cases the 
[Fe II/Mg II] ratio is consistent with meteoritic abundances). We also note the similarity 
between the velocity shift between the z=1.2234 system and the GRB 020813 host galaxy 
and the shift of the z=2.296 and z=2.298 systems with respect to the GRB 021004 host 
galaxy. Although the probabihty for a chance occurrence of the z=1.2234 system is not 
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as low as in the cases of the two shifted systems in GRB 021004, and that the z=1.2234 
ionization status is probably much lower than that of the z=2.296 and z=2.298 systems, 
this similarity might suggest a common scenario for the two GRBs. 

5. Conclusions 

One of the straightforward results of our UVES high resolution observations of two 
GRB afterglows is that the ISM of the host galaxies is complex, and many components 
contribute to each main absorption systems. These components span a total velocity range 
of up to about 3000 km/s. Several narrow components are resolved down to a width of a 
few tens of km/s. The UVES wide band coverage allowed us to investigate simultaneously 
both high ionization lines such as C IV and Si IV, and low ionization lines such as Mg II 
and Fe II in GRB 021004. This allowed us to constrain the ionization parameter of the 
gas of the differerent absorption systems. Combined with photoionization results obtained 
with CLOUDY the ionization parameters appear to lie in a relatively small range, with no 
clear trend with the system velocity. This can be interpreted as density fluctuations on top 
of a regular wind density profile. The [Mg II/Fe II] ratio of ^ 0.4 found for systems 
z=2.328_A, z=2.328_B, and z=2.298 is consistent with what is expected from meteoritic 
abundances. Assuming these abundances, the measured Si IV column density for system 
z=2.328_A implies a lower limit of ^5 x 10^^ cm~^ to the system total hydrogen column 
density. 

Indeed, our study shows that rapid reaction to the GRB triggers, high resolution 
and wide spectral coverage are the key ingredients to study GRB host galaxies. Today 
observations of this kind are still difficult and rather episodic. However, they should 
become routine after the launch of the Swift satellite, also thanks to the development of a 
dedicated Rapid Response observing Mode for the VLT telescopes. This mode will make 
possible automatic follow-up of GRBs (or other transient events) with response times in the 
lOm-lhr range, therefore helping in gathering spectra of unprecedented quality of medium 
to high redshift GRB host galaxies. 

We acknowledge support from contract ASI/I/R/390/02 and MIUR grant Cofin-2003- 
41. We thank Fabrizio Nicastro and Emanuele Giallongo for useful! discussions, Sandra 
Savaglio for her early work on this program, and an anonymous referee for comments that 
helped to improve the presentation. 
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